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0.1 Abstract 



Relative phase is treated as a physical quantity for two mode systems in quantum 
atom optics, adapting the Pegg-Barnett treatment of quantum optical phase to 
define a linear Hermitian relative phase operator via first introducing a com- 
plete orthonormal set of relative phase eigenstates. These states are contrasted 
with other so-called phase states. Other approaches to treating phase and pre- 
vious attempts to find a Hermitian phase operator are discussed. The relative 
phase eigenstate has maximal two mode entanglement, it is a fragmented state 
with its Bloch vector lying inside the Bloch sphere and is highly spin squeezed. 
The relative phase states are applied to describing interfcrometry experiments 
with Bose-Einstein condensates (BEC), particularly in the context of a proposed 
Hcisenberg limited interferometry experiment. For a relative phase eigenstate 
the fractional fluctuation in one spin operator component perpendicular to the 
Bloch vector is essentially only of order 1/N, so if such a highly spin squeezed 
state could be prepared it may be useful for Heisenberg limited interferometry. 
An approach for preparing a BEC in a state close to a relative phase state is 
suggested, based on adiabatically changing parameters in the Joscphson Hamil- 
tonian starting from a suitable energy eigenstate in the Rabi regime. 
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1 Introduction 



Studies of phase dependent phenomena in Bose-Einstein condensates (BEC) 
are hindered because phase has at least three different meanings. A similar 
situation applies in quantum optics [1] . In a first approach, phase is regarded as 
a physical property of the system [2] and is represented via a linear Hermitian 
operator that applies for all states and for which phase is the (real) eigenvalue. 
In a second treatment, the state is represented in a phase space [3J, [I] by a 
quasi-distribution function, and a complex phase is used to specify points in 
this space. In a third method, an operational approach [5], [6] emphasises an 
aparatus involved in measurements on the system and phase then refers to a 
feature of this measurement aparatus, such as the phase of a classical oscillator 
field interacting with the system. The dependence of system behaviour on phase 
in these three approaches are different in general, making comparison difficult 
because the meaning of "phase" is not the same. In a simple example presented 
in [I] , it is shown that no quantum state leads to a uniform phase dependence for 
all three different meanings for phase. It is important therefore to recognise that 
phase is different in the three approaches. Choosing which approach to use is 
somewhat a matter of personal preference, but from a fundamental point of view 
treating phase in the same way as other physical quantities would be preferable - 
if possible. The operational phase approach has the disadvantage of not linking 
phase to any intrinsic property of the system, and is dependent on the choice of 
aparatus. This generally involves some sort of homodyne system, but various 
combinations of beam splitters, phase shifters, vacuum input ports, detectors 
etc can also be involved. The phase space approach involves a complex phase 
which cannot be a measured value for any physical quantity and is dependent 
on the choice of distribution function used to describe the state. This could be 
Wigner W, Glauber-Sudarshan P, Husimi Q, .or other distributions. In fact the 
complex phase is often the eigenvalue of a non-Hermitian annihilation operator, 
which cannot represent a physical quantity. In neither of these two approaches is 
there any unique or compelling choice for defining phase. Thus the introduction 
of phase via eigenvalues of a linear Hermitian operator associated with the 
quantum system is the most objective approach pQ because it is not dependent 
on any particular way of specifying the state nor on any particular measurement 
system. This is not to claim that some ways of specifying the state are not more 
useful than others, nor is it intended to trivialise the difficult issue of measuring 
the probablity distribution for the measurable values of phase regarded as a 
physical property. Nor does it prove easy to find a suitable Hermitian operator 
to represent phase. However such an operator can now be defined both for 
quantum optical systems and Bose-Einstein condensates following the approach 
of Pegg and Barnett [5] that was originally applied to the quantum optics case 
(see [7] for a recent review). 

In the present paper we begin with a brief review of progress towards finding 
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a Hermitian phase operator, and then the Pegg-Barnett approach is adapted to 
define a relative phase operator for two mode systems via first introducing a 
complete orthonormal set of relative phase eigenstates. These states are con- 
trasted with other so-called phase states. Interesting properties of the relative 
phase eigenstates are then determined, these being entanglement, fragmentation 
and spin squeezing. In the final section applications of the relative phase states 
in describing BEC interferometry experiments are made based on treatment 
involving the Josephson Hamiltonian, and possibilities for preparing a BEC in 
a relative phase state are examined. Certain technical results needed for the 
main body of the paper are covered in an Appendix, which is available as online 
Supplementary Data. 



2 Hermitian phase operators 

2.1 Early attempts 

Early attempts to find a Hermitian phase operator for each mode by expressing 
the annihilation and creation operators a, at in terms of Hermitian number 
n a and phase operator <p a via a = exp(i(f) a ) (na) 1 / 2 , at = (n,,) 1 / 2 exp(— i<fi a ) 
with [<fi a , n a ] = — i [5], or by introducing two exponential operators [S], |10) via 
the annihilation and creation operators E + — (n a + l)~ 1,/2 a, E~ = at (n a + 

l) -1 / 2 = {E + y were unsuccessful. For the first, the commutation rule leads 
to a contradiction when matrix elements between number states are evaluated 
- (n — to) (n\(f) a \rn^ — —iS nm . For the second, the introduction of a Hermitian 

phase operator required being able to write E + = exp(i(j) a ) and E~ — exp(— i<f> a ) 
in the form of unitary operators also failed, for although E + E~ = 1 we have 
E~E + = 1 — |0) (0| - rather than 1 as is required. Note that Hermitian cosine 
and sine operators C, S can be introduced via E + = C + iS, but again these 
are not trigonometric functions of a Hermitian phase operator. However, the 
approach of Pegg and Barnett [5] via the introduction of phase eigenstates and 
then the Hermitian phase operator was successful. This approach does require 
introducing a cut-off on boson numbers for each mode, but this can be justified 
mathematically in terms of a sequence of Hilbert spaces [11] and shown not to 
affect physical predictions for finite energy fields. 

2.2 Relative phase operator and eigenstates 

In this section we introduce relative phase eigenstates and a Hermitian relative 
phase operator for the case of a two mode single component BEC with mode 
annihilation operators a, b and spatial mode functions <p a {r), <fib( r ) using a 
modification of the Pegg-Barnett approach for single modes. From the Fock 
state orthonormal basis states \n a ), \n ) involving n a , n bosons in the modes a 
complete orthonormal set of relative phase eigenstates \6 P ) for the N = n a + nb 
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boson system are then defined via 



N/2 

|0 P H 7 _= eMik0p)\N/2-k) a \N/2 + k) b (1) 

where P = p{2tt/{N + 1)), p = -iV/2, -iV/2 + 1, .., +7V/2 is a quasi-continuum 
of A^ + 1 equispaced phase eigenvalues. The Hermitian relative phase operator 
is then defined as 

e = 5>pl*p><*pl (2) 

v 

This approach has also been applied previously in [TJ], [TJ]. Also an un- 
normalised version of Eq.(p} with phase angle —9 P is introduced in [14] (see Eqs. 
A3 and A4). The relative phase operator defined in ([2]) depends on the choice of 
modes and the total boson number N, and its commutation law with the relative 

number operator frn = Wb-a^a) is [6, 5n\ = iJ2 K) (0 q \ (-l) p ~ q S^rm 

P^q 1 p q " 

rather than just i. The Fock state \N/2 — k) \N/2 + fc). is an eigenstatc of the 

relative number operator with eigenvalue k. Note that the present approach 

for a two mode system defines a relative phase eigenstate and relative phase 

operator, rather than phase eigenstates and operator for each mode. However, 

the relative phase operator can be defined without requiring a cut-off on boson 

numbers since there is an automatic restriction for k to lie between —N/2 and 

+N/2. The definition can be extended to apply to mixed state boson systems 

with a range of A" via Qt = Y^n®(N)Hn using projectors Iljv onto N boson 

JV 

states. Note that essentially the same states can also be defined for quantum 
optical systems, there the bosons are massless photons. 

An approach closer to the original Pegg and Barnett method would be to 
define phase operators for each mode, and then the relative phase operator would 
be the difference between the separate phase operators for the two modes, and 
this method is used in Ref. [15] . This approach requires introducing a cut-off on 
boson numbers for each mode and special techniques are needed to restrict the 
phase difference to a 2tt rather than Air interval. There are differences between 
this approach and that adopted here and in [TJ] , [TJ] , which are discussed in [TB] , 
[17] . The approach presented here provides a more direct focus on relative phase 
as a basic physical property and enables the relative phase to automatically lie 
in a 2tt interval. 

A similar approach to that here can also be used to define relative phase 
eigenstates and a Hermitian relative phase operator for a two component BEC 
where each (hyperfine) component is associated with a single spatial mode 
function. This situation is again a two mode system and similar Fock states 
to \N/2 — k) a \N/2 + k) b act as an othonormal basis, though now \n a ) has n a 
bosons in a spatial mode </> a (r) associated with internal (hyperfine) state a. 
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2.3 Pure states and quantum superpositions 

Any pure quantum state |<fr) for the N boson system can be expanded in terms 
of the relative phase states as 



N/2 

!<*>> = E A(e p )\e p ) 

p=-N/2 



(3) 



and the amplitudes A(6 p ) determine the probability P(9 p ) for measuring the 
relative phase 9 P via the standard expression 

P(9 P ) = \A(9 P )\ 2 (4) 

The same state can also be expanded in terms of the relative number states as 



N/2 

|$)= E b k \N/2-k) a \N/2 + k) b 

k=-N/2 



(5) 



with expansion coefficients b k . It is then easy to see that the expansion coef- 
ficients in terms of relative phase states and the Fock states are related via a 
Fourier transform. 



A ( P^ = /J-tt E exp(-ik6 p )b k b k = — L— . . 

V +1 k V N + 1 p 



'£exp{+ikO p )A(0 p ) (6) 



For the relative phase state itself the expansion coefficients are b k = exp(+ik9 p ) / y/N + 1. 
The generalisation for mixed states is straightforward. 

As an example of a quantum superposition we consider the state 



N 



N 
~2 



N N\. 



(7) 



which is the so-called NOON state, being a superposition of states | N, 0) and 
| 0, iV). It is also referred to as a Schrodinger cat state, and is an example 
of an entangled state. In the first term there are N bosons in mode 4>l and 
in mode 4>r and for the second the reverse applies. For this state b k = 

(5 k ,. 



-N/2 



$k,+N/2)/V2 and hence A(0 P ) = y/2/ (N + 1) cos(^# p ), which gives 
an oscillatory probability distribution for the relative phase with probabilities 
changing from to 2/(N + 1) for neighboring phase angles. Such oscillations 
would be hard to detect. On the other hand the different NOON state 



N 



N 



(8) 



with n <C N and A(0 P ) = y / 2/(N + 1) cos(n9 p ) would have a central peak in 
the phase probability for 9 p = and the first zero at P — ±n/2n 7 which would 
correspond to a relative narrow phase probability distribution with A9 P oc 1/n, 
if n is large enough. 



6 



2.4 Other phase dependent states 

Note that other authors [TB] have defined a set of states for a two mode BEC 
that depend on phase variables 9, \ via the expression 

\9, X ) = ^=( cos ^ ' expH^at + sin0 ■ cxp(+i± x fi) N \N) a \0) b (9) 

which are also referred to as phase states. For the case where 9 = ir/4 such 
states have been used to define a phase x hi BEC interferometry experiments 
[19] . and x i s measured in terms of the evolution time for a condensate in a 
double well trap when inter- well tunneling dominates over collisional effects. In 
this case phase is essentially the evolution time, which is an operational variable 
directly associated with the specific measurement process. However, states such 
as \9,x) are actually binomial states and correspond to all bosons being in the 
same single particle state cos9 ■ exp(—i^x) ( t ) a(' r ) + sin6* • exp(+i^x)4>b( r ) ■ They 
are also referred to as coherent states. Expanding the coherent states 

\0,X) =Ebk(e,X)\N/2-k) a \N/2 + k) b (10) 
k 

as a superposition of the basis states \N/2 — k) a \N/2 + k) b , the expansion co- 
efficients are bk(9,x) = ^N/2~k ' (cos 2 ~ fe ■ (sin9) N ^ 2+k ■ exp(ikx) involving 
binomial coefficients C^ /2 _ k = N\/((N/2- k)l(N/2 + k)\). The binomial states 
are physically important since they describe an unfragmented BEC [201 . How- 
ever they are not a complete orthonormal basis set for the two mode BEC. For 
example there is no choice of 9, x that gives the fragmented state |iV/2) a \N/2) b 
which has an occupancy for each of the two natural orbitals (see below) of 
N/2. For the coherent state with equal probabilities of finding a boson in each 
mode \9 = 7r/4,x) the expansion © in terms of relative phase states \9 p ) gives 
a relative phase probability 

for ./V large. This is a Gaussian distribution centred around 9 p = x with a narrow 
width of A9 P cx 1/y/N. For this coherent state the relative phase distribution 
corresponds to the standard guantum limit. Note that for large N this coherent 
state almost has a well-defined relative phase %, which may explain why it is 
sometimes regarded as being a state with a definite relative phase. However, 
they are not eigenstates of any relative phase operator. 

Other authors [5T] consider eigenstates \E(9)) of a phase dependent quadra- 
ture operator for each mode E(9) = i(aexp(—i9) — a* exp(—i9)) with eigenvalue 
E(9), and probability distributions given as | (E(0)\$) | 2 for finding the quadra- 
ture field to have an amplitude E(9) considered as a function of phase variable 
9 determine a phase distribution for a quantum state |$) without introducing 
phase as an eigenvalue of a Hermitian operator. The expansion of the quadra- 
ture eigenstates in terms of Fock states \n a ) a involves Hermite polynomials and 



7 



Gaussian functions of E(9). However, the states \E(6)) are non-orthogonal so 
the probability concept is doubtful. This treatment of phase is really an example 
of the phase space approach. 



3 Properties of relative phase eigenstates 

The relative phase eigenstate has several interesting properties. These include 
entanglement, fragmentation and spin squeezing. We deal with each in turn. 



3.1 Mode entanglement 

Firstly, it is a state with maximal mode entanglement for the a, b sub-systems, 
so is of interest in quantum information The entropy of entanglement is one of 
the standard measures of entanglement |22| and is given by the von Neumann 
entropy for the reduced density operator for either of the subsystems a or b. 
Thus for the system in pure state |$) the entropy of entanglement is 

S(f>a) = -k B Tr(p a \ogp a ) = ~k B Tr(p b \ogp b ) = S(pb) 

p a =Tr b (|$)($|) fc=Tr„(|3)<$|) (12) 

For the relative phase eigenstate it is straightforward to show that the entropy 
of entanglement is given by 

S(p a ) = k B \og(N + l) =S(fib) (13) 

which is very large. The general case of maximal mode entanglement occurs 
when the amplitudes b k in ([5J satisfy \b k \ = 1/y/N + 1 [23], and the relative 
phase eigenstate is a particular case. 



3.2 Fragmentation 

Secondly, it is a fragmented state |20j , since there are two natural orbitals with 
macroscopic occupancy. For large N the first order quantum correlation function 

G^\r, r') = /^(r^fr')^ (where ^(r), ^(r) are the usual field operators, 
\&(r) = a</> Q (r) + 606 (r)) is given by (see Appendix) 

G«(r,r') = ^;(r)0 a (O + ^(r)^(r')) 

+ ^exp(i0 p )(^(r)0 6 (r')) + exp(-z0 p )(0*(r)0 a (r')) (14) 

using the result that for large N the sum J (f (f + 1) - k{k ± 1)/((JV+1)) 
is approximately §iV. The natural orbitals are the eigenfunctions of the first or- 
der quantum correlation function, and are given by x±( r ) = ( ex P(^p/2)^a(r) ± 
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exp(— i9 p /2)4>l(r))/V2. The eigenvalues (which give the occupancies) are (| ± 
| )N. Fragmented states cannot be described by a single Gross-Pitaevskii equa- 
tion, generalised mean field theories are involved [53], [25] involving coupled 
generalised Gross-Pitaevskii equations. 



3.3 Spin squeezing 

Thirdly, the relative phase eigenstate is a spin squeezed state |26] in which one 
component of the spin angular momentum has a Heisenberg limited fluctuation. 
The Schwinger spin operators are defined by S x — (b^a + a%)/2, S y — (b^a — 
a*V)/2i, S z = (b^b — a^a)/2, and the Bloch vector is defined via its components 
(s x \ , (^Sy^ an d [27 , often in units of N. For large N the Bloch vector 

is determined to be (see Appendix) 

(s x ) = ^Ncos6 p (s y ) = ~Nsm9 p (d z ) = 

(15) 

using the result that for large N the sum £ fc \/(f (f + x ) ~ k(k ± 1)/((N + 1)) 
is approximately ^N = 0.3926iV. This vector is in the equatorial plane with 
azimuthal angle <f> = 2tt — 9 p , and is inside the Bloch sphere of radius N/2 - 
another indicator of fragmentation. For any unfragmented state is a coherent 
state \6, x) and the Bloch vector always lies on the Bloch sphere, the orientation 
being given by polar angle 7r — 29 and azimuthal angle 27r — \. Thus any state 
for which the Bloch vector lies inside the Bloch sphere must be a fragmented 
state, and the relative phase eigenstate is such a case. 

Spin operators along (J z ) and perpendicular (J x , J y ) to the Bloch vector 
may be defined by 

Jx = S z J y = S x sin 9 P + S y cos 9 p J z — S x cos 9 V — S y sin 9 p (16) 
and in terms of the new spin operators 

[j x )=0 (j v ) = (£\ = JiV a 0.392JV (17) 



The covariance matrix (see Appendix) which describes the quantum fluctua- 
tions for the spin operator components can be shown to be diagonal for the new 
spin operators J x , J y , J z . For large N the fluctuations in the new Bloch vector 
components for the relative phase eigenstate are found to be (see Appendix) 



SJ X fa y/T/UN « 0.289^ SJ y » J I + - hxN. SJ Z « y/(l/6 - tt 2 /64)N w 0.1127V 

V 8 4 

(18) 

where Sh 2 = ((fi - (&)) 2 ). As | (j z ) |/2 « 0.196N we see that for all N > 4 



the product SJ X ■ SJ y > 0.198AT, which is greater than \ (Jz) \/2 consistent with 
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the Heisenberg uncertainty principle. However, although J x is not squeezed, the 
other perpendicular component J v is highly squeezed, with a fractional fluctu- 
ation 8J y / (j*^ essentially of order 1/N due to the denominator (jzj- The 

numerator SJ y is a very slowly increasing function of N - for N changing from 
10 8 to 10 10 it only changes from 2.17 to 2. 42. The relative phase state could be of 
interest in Heisenberg limited interferometry |28j . By contrast, the fluctuations 
in the Bloch vector components for the coherent state \0,x) are SJ X ~ vN, 
5J y w y/~N and SJ Z w 0. Here the fluctuations are equal for the two compo- 
nents perpendicular to the Bloch vector, so there is no squeezing. Furthermore, 

the fractional fluctuation SJ XtV / (jz) is only of order 1/y/N, corresponding to 
the standard quantum limit and not to the Heisenberg limit, as is the case for 
the relative phase eigenstate. 



4 Applications of relative phase eigenstates 

The relative phase eigenstate provides a useful theoretical concept for describing 
interferometry experiments based on BEC. This type of application is discussed 
in this Section, firstly in general terms and then for a specific BEC interferometry 
proposal. However, before treating these applications the question of whether 
the relative phase eigenstate can be prepared via some sort of dynamical process 
will be examined. 



4.1 Creating relative phase eigenstates ? 

The energy and energy fluctuation associated with the relative phase eigenstate 
are quite large. The typical two mode system such as bosons in a double well 
potential is described via the Josephson Hamiltonian 

H = -JS X + SS Z + US 2 Z (19) 

where J is the inter-well tunneling parameter, S describes asymmetry of the 
two wells and U is the collision parameter. It is easy to sec that the relative 
phase state ((l) is not an energy eigenstate. The non-zero matrix elements of the 
Josephson Hamiltonian between the basis states \N/2, k) = \N/2 — k) \N/2 + k) h 
are 

H k , k =5k + Uk 2 



Hk.k+i = ~\j y(y + 1) - (k + l)fc 



tf fe , fe -l=-^y(y + l)-(fc-l)fc (20) 
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and for the relative phase state to be an energy eigenstate with energy E requires 



E = 5k + Uk 2 - y (y + 1) - (k + l)kexp(+i6 p ) 



J N.N 

-2Vy ( y + 1) 



(k — l)kexp(—i6 p ) 



(21) 



for all k, which is not possible. The mean energy (if ) is 



U ) = U—N 2 - Jcos9 v -N 
12 p 8 



(22) 



so in the Rabi regime |20j where J 3> UN the mean energy is approximately 
— JcosOp^N, whilst in the Fock regime [20] where U 3> JN it is essentially 



U^-N 2 

12 iv 



The variance in the energy SH = ( (i? — (H))) is given by 



£g2 ^ N 2 X ^ UN j s y 



l 

180 

COS ft 





384 w p 





" 




UN ' 





X 


J 


1 

12 - 
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. (23) 

correct to 0(N 2 ) (see Appendix for details). This is a quadratic form in the 
quantities UN, J and (5. That this form is positive definite can be shown 
by determining the eigen values Ai(0 p ), \2(0 P ) and Aa(# p ) of the 3x3 matrix 
in Eq. (|23p . and explict formulae are given in the Appendix. As expected the 
eigenvalues are all real and positive for all relative phase 9 P (see Figure A in Ap- 
pendix) . It is of some interest to consider cases where the Josephson parameters 
are related via 



(24) 



where [ X% a Xi a Xj, a ] are the orthonormal column eigenvectors associ- 
ated with the eigenvalues Ai(# p ), \i(9p) and \j,{9p) and K is arbitrary. In this 
case an expression for the relative energy fluctuation can be obtained as 



UN ' 




Xla 


J 


= K 


Xi a 


5 




X-ia 



H 



(Xia-p; - X- 



2a COS« P g / 



(25) 



for the eigenvalues \i(9 p ), A 2 (0 p ). For A 3 (6> p ) we have \H) = 0, so the relative 
fluctuation is undefined. In Figures 1 and 2 the relative energy fluctuations are 



shown for Ai( 



A, 



respectively. 
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Figure 1. Relative energy fluctuation for Josephson parameters in Ai(0 p ) case. 
Figure 2. Relative energy fluctuation for Josephson parameters in \2(0 P ) case. 

Clearly for the choice of Josephson parameters in the Xi(8 p ) case the rel- 
ative energy fluctuations are very large, 0(200%). However, for the choice of 
Josephson parameters in the \2{& P ) case the relative energy fluctuations are 
fairly small, 0(6%). If the Josephson parameters are chosen as in the latter 
case, then an adiabatic process starting with parameters as for some initial 9 p o 
and changing them to those for 9 p in accordance with Eq. ([24|) with a = 2 could 
prepare a state close to the required relative phase state. For example, with 
QpQ = f we have UN = 0, J = K and S — since X12 = 0, X22 = 1 and 
X32 = 0. A suitable initial state within the Rabi regime where J ^> UN, S = 
might be used. In particular, the state discussed below corresponding to that 
created at the end of the first stage in the proposed Heisenberg limited interfer- 
ometry experiment has a quite well defined relative phase 9 p q — (see Figure 4 
below) and might be suitable. The required quantities X12 and X22 that define 
the way UN, J would be adiabatically changed to reach any required 9 P are 
shown in Figure 3 (formulae are also given in the Appendix). 6 would remain 
equal to zero. 



Figure 3. Eigenvector for \2{0 P ). X\2 (blue curve) and X22 (red curve). 

A 32 = 0. 
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No actual experiment for preparing a BEC in a relative phase eigenstate 
or for directly measuring the relative phase probability have yet been carried 
out, so the relative phase probability distribution may need to be inferred from 
other measurements rather than directly measured. However, similar remarks 
may be made about position eigenstates for individual particles - where only 
states with relatively localised positions can be prepared and where position 
probabilty results are inferred from experiments involving scattering of weak 
probe beams, so this need not preclude the relative phase operator and its 
eigenstates being useful concepts in quantum atom optics. 

4.2 Interferometry experminents an quantum correlation 
functions 

The quantum correlation function G^(r, r') with r = r' is of particular interest 
as it determines the probability distribution for boson position measurements 
[29] . [30] and hence is useful in describing the interference fringes that can 
occur in BEC interferometry experiments. For a general state ([3]) the first order 
quantum correlation function can also be expressed in terms of the amplitudes 
A(9 p ) for the relative phase eigenstates. These amplitudes appear via three 
autocorrelation functions. We have 

G«(r,r') 




(26) 

where the autocorrelation functions of the amplitudes A(9 P ) are defined as 

N 

C Q (9 r ) = ± A({9 r + 9 q } mod27T ) A(9 q )* 
?=-¥ 

N 

C±i(6 r )= J A({0 r + q } moA2 „) A(6 q y exp(±iO q ) (27) 

Hence we see that for a state with a relative narrow relative phase distribu- 
tion around a particular phase 9q (A(9 p ) ~ ^e p ,e ) the autocorrelation func- 
tion Co(9 r ) will be peaked around 9 r = 0, whilst the autocorrelation functions 
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C±i(9 r ) will be peaked around 9 r = =F#o- This means that the first two terms in 
G( 1 )(r, r') have no dependence on O , whereas the last two terms have essentially 
a sinusoidal variation with 9q since the sum over 9 r will be dominated by terms 
&r =f T@o- If the particular central phase 9q is changed during an experiment 
then the boson position probability will change - hence a fringe pattern would 
be observed. Note that the observation of the fringe depends on the overlap 
of the mode functions being sufficiently large. On the other, for a state with 
a relatively wide relative phase distribution the auto correlation functions will 
be significant for a wide range of 9 r so the sum over 9 r will be no longer be 
dominated by terms 9 r = =F#o an d the fringe pattern would be washed out. 

4.3 Heisenberg limited BEC interferometry experiment 

The relative phase eigenstate is a valuable theoretical concept for describing the 
behaviour in BEC interferometry experiments. For example, in the proposed 
experiment by Dunningham and Burnett [3T] for Heisenberg limited interferom- 
etry in two mode BEC, the collapse and revival of interference fringes can be 
discussed in terms of collapses and revivals of the time dependent probability 
distribution for the relative phase. Collapse and revival effects in BEC were 
described earlier by Wright et al |32j . 

The Dunningham and Burnett experiment treats the two mode double well 
BEC system via the Josephson Hamiltonian. The experiment has two stages. 
In the first stage the system starts with equal numbers of bosons in each well, so 
the quantum state is |$(0)) = \N/2) a \N/2) b , so & fc (0) = S k ,o- From Eqs.Q and 
(O it is easy to see that the relative phase probability distribution is uniform. 
With evolution dominated by the tunneling term the state evolves for a time 
Xi = nh/2J (or when (f> = Jt/h = tt/2) The methods of angular momentum 
theory can be used to determine the dynamics, since the evolution operator 
U(t) — exp(iS x Jt/h) is just a rotation operator. We find that 

b k {T x ) = expOfc-) (jv/2)!(2jv/2) 2J-i; ^ C P+k C P 1 (28) 

from which we can calculate the relative phase probability distribution via 
Eqs.(j4]) and (|6]). This is shown in Figure 4 for the case N = 80 and we see 
that the system has a well-defined relative phase of approximately zero. This 
stage of the experiment involves creating a state with a rather well-defined rel- 
ative phase. Experimentally the time T\ is determined by observing the time it 
takes for the fringe pattern to become sharpest. 



Figure 4. Relative phase probability at end of the first stage. Well defined 
relative phase seen. Parameters are given in text. 
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The second stage involves evolution for a further time T dominated by the 
collision term, or the collision term plus the asymmetry term. In this case we 
find that 



ST 

bkiTx + T) = exp(-ifc— ) exp(- 
n 



_ zfc 2^ )MTi) 



We first consider the situation when there is no asymmetry S = 0. A character- 
istic time scale for the collision dominated evolution is T 2 = Trfo/2U (or when 
£ = Ut/h = 7r/2). However, due to the exp(— ik 2 ^-) factor there is a dephas- 
ing effect, causing the relative phase probability amplitudes A(9 p ,Ti + T) to 
become significant over a wide range of 9 p . This causes a collapse in the previ- 
ously well defined interference fringe pattern. The time scale for this to happen 
is that required for the fastest pairs of contributions (k = —N/2,— N/2 + 1 or 
k = +N/2, +N/2 — 1) to get out of phase by ~ it. Thus the collapse time is 
given by T c = irh/NU = T 2 /N (or when f = ir/2N), which is 0(1/N) times 
shorter than T 2 . This collapse effect is shown in Figures 5 and 6 for N = 80. 
For £ = 0.0017T the relative phase distribution is starting to spread out and is 
essentially uniform when £ = 0.017T. 



Figure 5. Relative phase probability just after the end of the first stage. 
Dephasing effects starting to be seen. Parameters are given in text. 

Figure 6. Relative phase probability somewhat after the end of the first stage. 
Complete dephasing effects seen. Parameters are given in text. 
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However, the factors exp(— ik 2 ^^) do eventually get back into phase. If 
£ = Ut/His a multiple of ir/2 then all the phase factors have a modulus of unity, 
irrespective of k. Hence a revival of the relative phase probability distribution 
to the sharply defined distribution that occured at the end of the first stage 
will take place. The revival time scale is thus given by T rev = T 2 — ttH/2U (or 
when £ = Ut/H = 7r/2). This is shown in Figure 7 for N — 80. In Figure 8 the 
time is slightly longer than T rev and the relative phase distribution is starting 
to collapse again. 

Figure 7. Relative phase probability at the end of the second stage. Revival of 
well-defined phase seen. Parameters are given in text. 

Figure 8. Relative phase probability just after the end of the second stage. 
Beginning of collapse of well-defined phase seen. Parameters are given in text. 

We now consider the effect of asymmetry. It is easy to see that at time 
T\ + T 2 the relative phase amplitude for non zero 8 is given by 



so is of the same form as when there is no asymmetry, but with the relative phase 
shifted by 5T 2 /h = (5/U)n/2. This effect is shown in Figure 9 for N = 80. The 
shift in the fringe pattern would be observable if 8 is a reasonable fraction of U. 



A(9 P , Ti + T 2 ) = A{8 P + ST 2 /H, T x + T 2 ) s=0 
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Figure 9. Relative phase probability at the end of the second stage. 
Asymmetry present. Well-defined but shifted phase seen. Parameters are 

given in text. 

Evolution during the second stage of the experiment is allowed to occur for 
a time Ti corresponding to the revival time, and with zero asymmetry present. 
The revival time could be determined experimentally by observing when the 
sharp fringes obtained at the end of the first stage are restored again. The 
accuracy in determining the revival tiime is given by the collapse time, so the 
fractional error in the revival time T rev is of order 1/N. If the second stage 
is run again with asymmetry present the fringe pattern at the revival time is 
shifted by 8T rev /h. If this phase shift is measured with perfect accuracy, then 
the fractional error in measuring S is the same as that for T revi and hence is of 
order 1/N. This represents a Heisenberg limited interferometry measurement 
of the asymmetry scaling as the inverse of the total number of bosons. 

5 Summary 

This paper presents an approach to treating phase in quantum atom optics in 
which phase is regarded as a physical quantity for the system and treated theo- 
retically as a linear Hermitian operator, similar to the Pegg-Barnett treatment 
of phase in quantum optics. Other approaches to treating phase are discussed 
and a brief review outlines previous attempts to find a Hermitian phase opera- 
tor. The Pegg-Barnett approach is adapted to define a relative phase operator 
for two mode systems via first introducung a complete orthonormal set of rel- 
ative phase eigenstates. These states are contrasted with other so-called phase 
states. The entanglement, fragmentation and spin squeezing properties of the 
relative phase states are set out. The relative phase state has maximal two 
mode entanglement, it is a fragmented state with its Bloch vector lying inside 
the Bloch sphere and is highly spin squeezed. In the final section applications 
of the relative phase states in describing BEC interferometry experiments are 
made, both in general and in the context of a proposed Heisenberg limited in- 
terferometry experiment. Interferometry experiments essentially measure the 
autocorrelation functions for the relative phase amplitudes. The possibility for 
preparing a BEC in a relative phase state is examined, and an approach based on 
adiabatically changing parameters in the Josephson Hamiltonian is suggested. 
However, the difficulty is similar to preparing a particle system in a position 
eigenstate. In spite of this, the relative phase states are still a useful concept 
for describing experments in quantum atom optics. Finally, if such a highly 
spin squeezed state could be prepared it may be useful for Heisenberg limited 
interferometry in view of the fractional fluctuation in one of the spin operator 
components perpendicular to the Bloch vector being essentially only of order 
1/N. 
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7 Appendix 

In the appendix we summarize certain key results needed in the main body of 
the paper. 

7.1 Field operators 

For the two mode BEC the field operators are given in the two mode approxi- 
mation as 

$(r) = acf> a (r) +Hb{v) *(r)t =tf<p a (r) +V<fi(r) (29) 

and with the usual non-zero commutation rules for the mode operators [a, a 1 '] = 
[b, &t] = l,the non-zero commutation rule for the field operators is 

[$(r), ^(r'jt] = ^(r)^(r') + b (r)^(r') = 6 2 (r, r') (30) 

which is a restricted delta function for the space spanned by the two orthonormal 
mode functions. 



7.2 Spin operators and spin states 

Because of the two-mode approximation it is possible to treat the bosonic system 
using the methods of angular momentum theory. The system behaves like a 
macroscopic spin system with angular momentum quantum number j = ^ . 

In a two-mode theory it is convenient to introduce the Schwinger spin angular 
momentum operators defined by 

S x = ($a + a^b)/2 
S y = (&a-affi)/2i 

S z = (pb-a^a)/2 (31) 

The spin operators S a satisfy the standard commutation rules for angular mo- 
mentum operators 



= i£ a bcS c (a,b,c = x,y,z), (32) 



and the square of the angular momentum (S^) 2 can be related to the boson 
number operator. Thus: 

(|) 2 -E(^) 2 
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where 



N = (pb + a^a) 



(34) 



is the number operator. Clearly the angular momentum squared is a conserved 
quantity. Note that the spin operator S z is the same as the relative number 
operator SN. 

The N boson system behaves like a giant spin system in the two-mode ap- 
proximation. The basis states \N/2 — k) a \N/2 + k) b are simultaneous eigen- 
states of (S^) 2 and S z with eigenvalues + 1) and k respectively. To em- 

phasize the spin character of the basis states we can introduce the notation 



\N/2-k) a \N/2 + k) t 



(f+fc) 



m -*)!]* [(t 



2'* 



0) 



(35) 



Thus: 



(|) 2 



N 



k) = 



2 ' 



N N 
~ ~2^~2 



= k 



1) 



2 ' 



N 



(36) 
(37) 



Y is the spin angular momentum quantum number, and k is the 



Hence j 

spin magnetic quantum number, with (— -j < k < -y). Thus the boson number 
N and the quantity k that specifies the fragmentation of the BEC between the 
two modes have a physical interpretation in terms of angular momentum theory. 
Since boson numbers may be ~ 10 s the spin system is on a macroscopic scale. 

As in angular momentum theory we find it convenient to introduce spin 
up 5+ and spin down S— operators, which change the spin magnetic quantum 
numbers by ±1. We have 



S± 



N 



N N 

{-(- + i)-Hk±i)y. 

N N i 

{(- T k)(-±k + l)}-* 



S± — S x i iSy 



N 



N 

y 



,/c± i 



k± 1 



(38) 



The methods of angular momentum theory can be utilized by first writing 
the full Hamiltonian in terms of spin operators using equations (|2"5|) . (|3"Tj) , (|34p 
- noting that all terms involve equal numbers of creation and annihilation oper- 
ators, and its matrix elements calculated using angular momentum theory from 
(|3"T|) and (|3"5)) . The same applies in a simplification to the full Hamiltonian 
giving rise to the Josephson Hamiltonian. 
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7.3 Quantum correlation functions 

The first order quantum correlation function is defined as 

G (1) (r,r') = ^(r)**^')) (39) 
and for the pure state given by ([3]) this is found to be 

G«(r; r') = £ b k *b k |^ (r)> a (r') - kj + Mr)* Mr') (f + *) } 



+ E&fc*6fc-i|^WVa(r')y(y + fc) (f- fc + 1 )j ( 4 °) 
7.4 Bloch vector 

The components 5 a of the Bloch vector are given by averages of the spin oper- 
ators S a 

S a = (s a ) (a = x,y,z) (41) 



Often the Bloch vector components are scaled in units of N, but to avoid extra 
notation we will not do that here. 



Since S a is hermitian and (S a ) < ({S a ) 2 ) and using (|33|) we see that 



o<£s*<£((s a ) 2 ) = f(f + 1) = ^ N>>1 (42) 

a a 

showing that for all states the Bloch vector lies inside or on a Bloch sphere, 
whose radius is -y. 

For the quantum state given by (JSJ expressions for the Bloch vector compo- 
nents are 



s± = E^Tiy (f (f + 1) - m t i) 

5, = (S+ + S-)/2 S y = (5+ - S_)/2z 

S, = (43) 
k 

and these are related to the first order quantum correlation function via 
G (1) (r;r') 

f {MrrUr') + Mr)*Mr')} 
+S Z {-Mr)* Mr') + Mr)* Mr')} | (44) 
+S- {Mr)* Mr')} + S+ {Mr)*Mr')} 
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7.5 Covariance matrix for spin operators 

The covariance matrix C(S a ,S b ) for the spin operators S a is given by 



C(S a , S b ) = l - ((A5„ AS b ) + (AS b AS„)) (45) 
AS a = £ a -(s a \ (a, b — x,y, z) (46) 



where A,S a is a spin fluctuation operator. It is easy to see that the 3 x 3 co- 
variance matrix is real and symmetric and that C(S a ,S a ) gives the variance 

(^(AS a ) 2 ^ for S a . These are the square of the standard deviations or fluctua- 
tions. Such a matrix defines a positive quadratic form F(^ x ,S_y,^ z ). With real 
£ a we have 

a, b 

= \ (^E£° A ^X> A ^ + ^E& A ^X> A ^) 

= (S0S(O) > (47) 
for any state, where S(£) = ^(aA^ = S(^. Hence the three eigenvalues for 

a 

the covariance matrix will be real and positive. Linear combinations of the AS a 
involving a real orthogonal matrix will diagonalise the covariance matrix and 
the diagonal elements will give the variances for fluctuations in three orthogonal 
directions. These specify the principal quantum fluctuations. 
The covariance matrix can also be written as 

C(S a , S b ) = ^({s a S b + S b S a ) ) - (s a ) (s b ) (48) 

so it measures the difference between the average of half the anti- commutator 
of S a , S b and the product of the averages of the separate S a ,S b . 

Expressions for the covariance matrix elements for the spin operators in the 
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case of the pure state given by ([5]) are as follows. 
C 

*~^xx 

= C(S X , S x ) 



-JZ b *^b k \J y (y + 1) - Hk + j(j + 1) - (fc + l)(fc + 2) 

A; 
fc 

lE^ 6 



2 v 2 



AT N N N 

-(IT + !) - - + !) - ( fc - - 2 ) 



2 v 2 



/ YlK+ihy/m + l)-k(k + l) \ 
k 

V fc / 



(49) 



and 



C{S X , Sy) Cya; 



4/ 



2 v 2 



N N N N 

i W -(-S- + 1) - *(* + + !) - ( fc + + 2 ) 



2 v 2 



iV iV N N 

, w _(_ + 1) _ fc fc - lW-(- + 1) - (fc - l)(fc - 2) 



M/ ^(^ + i)-fc(fc + i) 



2 v 2 

2 



(50) 
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and 



= C(S X , S z ) = C z 



^K+MZk + 1) sj y ( y + 1) - + 1) 



+ ^K-M^ 1) \j y ( y + 1) - fc(fc - 1) 



lT,<>UihyJ y (y + 1) - Mfc + fc 
l^k-ihyj y (y + 1) - - l)5>fc 6 * fc 



and 



C(Sy,Sy) 



k+2 b k 



N ,N 



2 v 2 



fe 



AT, AT 



(51) 



(- + !)- k(k + lW -(- + 1) - (k + l)(fc + 2) 



2 v 2 



JE 6 *-A V f (y + !) - - 1) V y (y + 1) - (k l)(fc - 2) 



2 v 2 



1 

+ 4 



fe 

-YiK-ib^m +i)-fc(fc-i) 



(52) 
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and 



a 



C(S'j / , S Z ) — C Z y 

= p^6^ +1 6 fc (2fc + 1)JE Z + i)_fc(fc + i) 



1 x-^ /TV JV 

-^6^1^(2^ - 1) J -( T + 1) - - 1) 



(53) 



and finally 



Czz — C(S Z , S z ) 



k \ k ) 



(54) 



In terms of the new spin operators defined via the orthogonal transformation 
in (|16p . the new covariance matrix is given by 



C(J a ,J b ) = Y, M °c(.0p)C(S c , S d )M bd (dp) = Sab ( (A J ) 



(55) 



where 



[M{6 P )] = 



sin 6 
cos ( 



sin 9 P 



cos ( 



relates the new and original spin operators via 

Ja = ^~]Mac(9p) S c 



(56) 



(57) 



It turns out that the new covariance matrix is diagonal. The evaluation of the 
original covariance matrix involves the following sums for N large: 

(a) £ fc X /(f(f + l)-fc(fc±l)/(JV + 1) = z£ 



GO £ fc V(f (f + 1) " Hk ± 1)V (f (f + 1) " fa ± 1)(* ± 2)/(AT+l) = f- 
(c) E fc (f(f + l)-fc 2 )/(7V + l)=^ 
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(d) £ fc (2ft±lW(f (f + 1) - k(k ± l)/(N + l) 4= 0. These sums are correct 



to 0(N). This gives the new covariance matrix to 0(N) as 
C(J a ,Jb) = 










6 64 1 Jv 



(58) 



showing that correct to O(N) the variances for J. x , J z are large but that for J y 
is zero. 

To determine the new covariance matrix more accurately we begin again 
with the new spin operators still given by (|57[) and work out the expressions for 
C( J a , Jb) to O(N ) when the term is zero correct to 0(N). As before we find 
the new covariance matrix is diagonal but now given by 



C(J a , Jb) 



(T2) N2 






I In AT 



\6 64 J 



N 2 



The new variance for J y is now found to be non-zero and given by ~ 
This requires the following sum: 



(59) 



■In AT. 



(e) (- Efe V(f (f + 1) - Hk ± 1) J(f (f + 1) - (fc ± 1)(* ± 2)/2(JV+l)+ 



(Kill 

l^k\ 2 V 2 



fc(*±l))/2(2V+l)) 



ilniV. 



7.6 Energy fluctuations 

The Josephson Hamiltonian is given by 

H = US 2 - JS X + 8S Z (60) 
and it is straightforward to show that then variance in the energy is given by 



SH 2 = [ U -J S ] 





U 


X 


-J 




s 







{S X ,S X } {S X ,S Z } 
{S Z ,S 2 } {S Z ,S X } {S Z1 S Z } 

( 61 ) 

which involves the covariances of S 2 , S x and S z . Evaluating the covariances 
using certain sums given in the previous section plus 
(f) £ fe k 2 /(N + 1 ) = ^ (g) E k kV(N + 1) = ^ 

(h) J2k H k + WCf (f + 1 ) -k(k + l)/(N + 1) = w we have correct to 
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0(N 2 ) 
SH 2 

= [u 



-J 5 



-J- N 4 

180 



= N 2 x [ UN J 5 ] 



i 

180 

384 COS ®P 



384 
'I _ Hi 
>6 64 



TV 3 cos ft 












" U 





X 


-J 


J-N 2 

12 JV 




5 



3k COs9 P 

(i-£) C os 2 e p 






_1_ 

12 



UN 
J 
5 
(62) 



which is a quadratic form in the quantities UN, J and S. That this form is 
positive definite can be shown by determining the eigen values Ai(ftp), A 2 (ftp) 
and X^iOp) of the 3x3 matrix. 

The corresponding eigenvalue equations are 



l 

180 

384 C0S 9 P 



384 cos @p 

(|- f^)cos : 





" 




Xi a 









X 


Xi a 


— A a (ftp) 


X^a 


1 

12 - 




X^a 







(63) 



where the eigenvectors are orthogonal and normalised to unity ^ Xi a Xip = 

i 

5 a p. The eigenvalues are easily obtained as 



Ai(ftp) = g ((a + b ) + V(a-b) 2 +4c 2 ) 



A 3 (ftp) 



((a + b) - v/(a-6) 2 +4c 2 ) 



(64) 



with a = y|q, b — (| — |j) cos 2 ftp, c = ^ cos# p and = ^. The eigenvalues 
are all real and positive, as may be seen in Figure 10 for the non-trivial Ai(ftp), 
A 2 (ftp). 



Figure 10. Eigenvalues Ai(ftp) (blue curve) and \2{0 p ) (red curve) 
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Eigenvalues 



0.015 



0.010 



0.005 



The corresponding eigenvectors are: 

X12 
X22 
X32 

X13 
X23 
X33 



V(a-A!) 2 + C 2 
1 

V(o - A 2 ) 2 + c 2= 

" " 


1 



(a- Ax) 


— c 

(a-A 2 ) 




The eigenvector for \i(9 p ) are shown in Figure 11. 



Theta 



(65) 



Figure 11. Eigenvector for Xi(6 p ). Xn (blue curve) and AT 2 i (red curve). 

*31 = 0. 



For the case where 



UN ' 




X\ a 


J 


= K 


X2a 


5 




X3 a 



(66) 



32 



EigenVectorOne 
1.0 h 



0.5 




where K is arbitrary, it is easy to see that in this case the energy variance and 
standard deviation are given by 



N 2 K 2 X 



5m) =NKJ\ a (6 p ) 



(67) 



The average energy in this case is 



II) = A ( UN— - J cos 6» - 
12 p 8 



//) = -V/v (X 1q ^-X 2q cos0 p ~ 



(68) 



so the relative energy fluctuation is 



y/X a (6 p ) 



| (-^la — -^2a COS #p f ) 



(69) 



This expression only applies to the Xi(0 p ), A 2 (6 ) p ) eigenvalues, since for X 3 (9 p ) 



we have = 0. 
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